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ABSTRACT 

We look for the expected signature of an accretion disc by examining the properties of 
the Ha emission line versus viewing angle in a sample of 22 superluminal (SL) quasars. 
The Doppler factor 5 , jet velocity 7, and viewing angle towards the jet 9, are derived 
from published radio and X-ray data. Most of the Ha spectra (fourteen) have been 
observed at the United Kingdom Infrared Telescope ( UKIRT ) and are reported here. 
About one fourth of the SL objects have weak or absent Ha emission lines, with small 
equivalent widths (EW). These have high optical polarization, radio core dominance, 
Doppler factor and most of them high apparent SL velocity and low viewing angles. 
Therefore these weak-EW objects almost certainly have relativistically beamed optical 
continua. The strong-EW objects also show a clear beaming effect, but a much weaker 
one, with line EW varying by only a factor three while radio core dominance varies 
by a factor of several hundred. The correlation of EW with 9 is quantitatively in good 
agreement with the prediction of a flat accretion disc with limb darkening. The weak- 
and strong-EW sources also show an anti-correlation of line velocity width with the 
various beaming indicators. Again, the correlation with the derived viewing angle 9 
shows a quantitative agreement with the effect expected for an axisymmetric structure 
with velocity dominated by rotation. The line emission cannot come from the surface 
of the disc, or the line beaming would cancel the continuum beaming. However it could 
come from an axisymmetric system of clouds co-rotating with the accretion disc. 

Key words: galaxies: active - quasar: emission lines - quasars: general - radio contin- 
uum: galaxies - accretion disc 



1 INTRODUCTION 

Attempts to incorporate all AGN (active galactic nuclei) 
activity within some framework ("Unified Schemes") have 
focused on the importance of the orientation of the ob- 
server with respect to the source, in three separate ways. 

(i) The most dramatic effect is relativistic Doppler-beaming 
in a jet of radio-emitting plasma aligned close to the line of 
sight. Amplification of the apparent brightness of the source 
and apparent superluminal motions are then expected (Rees 
1966; Lind & Blandford 1985; Ghisellini et al. 1993; Urry 
& Padovani 1995). This effect explains the existence of 
"blazars" , objects dominated by a polarised and highly vari- 
able continuum, as well as one sided radio jets and other 
phenomena. This is the most well established orientation ef- 
fect, but only applies to the minority of radio-loud objects. 

(ii) The second effect is that of absorption or obscuration, 
most probably by an equatorial disc or torus, invoked to 
account for the dissimilarities between type 1 and type 2 
Seyfert galaxies (Lawrence & Elvis 1982, Antonucci & Miller 



1985, Antonucci 1993). (iii) Finally a third orientation effect 
is associated with an accretion disc geometry: (a) a simple 
optically thick, geometrically thin disc viewed inclined will 
appear fainter than when it is viewed pole on. However rela- 
tivistic effects in the vicinity of the accreting black hole make 
edge-on discs brighter (Cunningham 1975). (b) If the emis- 
sion lines share the same geometry they will have a broad 
velocity width at low latitude, (c) The polarization of the 
disc radiation will increase with the inclination (e.g. Rees 
1975). 

Some authors have searched for evidence of this third 
class of orientation dependent effects, generally by using the 
ratio of core radio luminosity to extended radio luminosity 
(e.g. Orr & Browne 1982; Wills & Browne 1986; Browne & 
Murphy 1987; Jackson et al. 1989; Jackson & Browne 1991b; 
Corbin 1997). This ratio is an orientation indicator (see for 
instance Ghisellini et al. 1993, Wills & Brotherton 1995) and 
these studies have been consistent with the widespread pres- 
ence of optical beaming in radio loud quasars in general (as 
opposed to simply in blazars, e.g. Browne & Murphy 1987) 
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and a flattened geometry for the emission lines (e.g. Wills & 
Browne 1986, Corbin 1997). The amount of optical beaming 
seen is far less than seen in the radio regime, indicating that 
it might be due to to relativistic beamed continuum seen 
in low contrast against less anisotropic disc emission (e.g. 
Browne & Murphy 1987; Wills 1991; Wills et al. 1992). The 
idea of this paper is to use superluminal sources for which we 
can derive an actual viewing angle rather than an indirect 
orientation indicator, and so directly test for the quantita- 
tive effects of accretion disc surface brightness and emission 
line velocity projection (see also Wills & Brotherton 1995). 

In this paper, we study the Ha emission of a sample 
of 22 superluminal sources. This line was chosen because it 
is the strongest low ionization emission line (LIL, e.g. hy- 
drogen Balmer, Mg II etc.) for which, there are both obser- 
vational (e.g. Wills & Browne 1986) and theoretical (e.g. 
Collin- Souffrin 1987) suggestions that its velocity field is 
anisotropic, likely in a rotational disc structure. In addi- 
tion, Ha line is less affected by absorption or blended opti- 
cal Fe II multiplet emission than H/3 . Line width is used to 
study velocity projection, and equivalent width to indicate 
the brightness of the underlying continuum. The viewing 
angle is derived from published SL motions, radio and X- 
ray fluxes, and angular sizes, using the method of Ghisellini 
et al. (1993). For eight sources we use published Ha data, 
but for the majority we have obtained new data. Most su- 
perluminal sources are at moderately large redshifts, so this 
has required obtaining IR spectra, which we report in this 
paper. 

In Section 2 we summarise the technique used to de- 
rive viewing angle. In Section 3, we describe the SL sam- 
ple giving relevant data and the parameters derived by the 
Inverse Compton model. In Section 4, we present the new 
UKIRT observations, with detailed source notes in Section 
5. In Section 6 we present the statistical results which are 
discussed in Section 7 and summarised in Section 8. The 
cosmological model assumed in this paper is H = 75 km 
s _1 Mpc -1 and g o =0.5. 



2 METHOD OF DERIVING VIEWING ANGLE 

Derivation of the viewing angle relies on two effects - su- 
perluminal motion, and Doppler boosting. We follow the 
method of Ghisellini et al. (1993) which we now summarise 
briefly. 

Apparent Supcr-Luminal (SL) motions of individual ra- 
dio structures (measured with VLBI experiments, see for 
instance Zensus & Pearson 1987) have been observed in ex- 
tragalactic radio sources for more than 20 years. These mo- 
tions are usually modelled by one or more radiating "blobs" 
moving at relativistic velocity /3 (in units of c) away from 
a stationary "core" . If the viewing angle is 9, the apparent 
transverse velocity is 



and can exceed the true speed and even that of light 
if the direction of the motion is close to the line of sight. 
The observed SL velocity therefore gives a combination of 
the Lorentz factor 7, and jet viewing angle 9 . The radiating 
material is brightened by the Doppler beaming factor 8 , 
which also provides a combination of 7 and 9 (8 = [7(1 — 
P cost))] -1 ). If 8 can be separately estimated, then, following 
Appendix B of Ghisellini et al. (1993), 7 and 9 are given by 



/3 a 2 + S 2 + 1 
28 
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The value of 8 is estimated from observations of the flat 
spectrum radio core, assumed to be made up of multiple self- 
absorbed components, by requiring that the Inverse Comp- 
ton flux of the component peaking at the observed radio 
frequency is less than the observed X-ray flux (see equation 
1 of Ghisellini et al. ). This also requires the angular size of 
the source, from VLBI measurements, and an assumption 
for the underlying synchrotron spectral slope, taken to be 
a = 0.75. 

We summarise uncertainties in the estimate of 9. (i) The 
value of 8 estimated from the X-ray flux is a lower limit, and 
so the derived 9 is an upper limit. However the value is very 
insensitive to this effect - 8 oc F^T 0,22 . So if the true F x is 
actually lower by a factor 10, 100, or 1000, the value of 8 is 
larger by a factor 1.5, 2.3, and 3.5 respectively. This factor 
might increase by up to 30 % if the source is not spherical as 
assumed, but a continuous jet (see equation 2 of Ghisellini 
et al. ) (ii) Radio flux and angular size errors may also in- 
troduce uncertainties on 9, but in most cases these are not 
large, (iii) Another possibility is that the superluminal ve- 
locities represent a pattern velocity which is different than 
the bulk velocity (Lind & Blandford 1985; Vermeulen & Co- 
hen 1994). (iv) Finally the observed value of /3 a and thus 9 
depends on the assumed Hubble constant, but we take that 
to be understood (Freedman et al. 2001). If H is smaller 
than the assumed value, 9 is smaller in most cases. In sec- 
tion 6 we provide arguments that 8 cannot be wrong by a 
very large factor, but also illustrate the effect on derived 
values of 9 of assuming that 8 is wrong by a factor of a few. 



3 THE SAMPLE 

The sample of the SL sources was selected from the list of 
Vermeulen & Cohen, 1994 (VC94). The final sample used in 
this paper is the superset of those for which we were able 
to obtain good quality IR spectra at UKIRT , plus those 
with published (optical) data on Ha . This sample repre- 
sents more than the one third of the known SL extragalac- 
tic sources. It is listed in Table 1. Columns (1), (2), (4), 
(5), and (6) list names, redshift, the fastest measured ap- 
parent velocity, and the core-to-lobe flux ratio R at 5 GHz 
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as given in VC94. Column (3) divides the objects into three 
classes depending on radio core dominance and optical po- 
larization. Core dominated sources with low polarization are 
shown as "CD" and those with high polarization (P > 3%) 
as "CDP" . Lobe dominated sources are shown as "LD" . The 
limit between CD and LD was taken as log R =-0.02. 

Columns (7)-(10) give the observed radio core flux and 
angular size {6d = Vab where a and b are the major and 
minor axes in mas) at a given frequency, together with the 
reference. Finally the two last columns give the observed X- 
ray flux at 1 keV together with the reference. The radio and 
X-ray data were then used to estimate 8 , and thus 9, and 7, 
as explained in the section above. For 3C 395 we have not 
found any X-ray measurement published, and we have used 
optical data instead. Table 2 shows the result, with sources 
listed in order of increasing 9, which can be seen to correlate 
clearly with S but not with 7, as expected in the beaming 
model. 



4 UKIRT SPECTRA 

4.1 Observations and Data reduction 

From the VC94 sample only objects between RA 05^ and 
19 h and Dec. lower than 60° were observed, due respectively 
to the time of year of the observations and the declination 
limit of the telescope. Finally we selected sources at red- 
shifts that would place Ha well inside one of the J, H, or K 
atmospheric windows. This gave 14 sources which are listed 
in Table 3. 3C 286 (not known as superluminal) was also 
observed during an RA gap in the run and is likewise listed 
in Table 3. 

The spectra were obtained with the common-user 
Cooled Grating Spectrometer CGS4 at the f/36 focus of 
the 3.9-m UKIRT on Mauna Kea, on three consecutive pho- 
tometric nights starting on UT date 1994 March 27. The 
instrument was equipped at this time with a 62x58 InSb 
array detector at the focal plane of a 150-mm camera with 
a slit of width equal to 3 arcsec, stepped perpendicular to 
the dispersion. The detector was stepped in increments of 
one-third of a resolution element in order to properly sample 
the instrumental profile. Different gratings and orders used, 
depending on the redshift of the source, and exposure times, 
are given in Table 3. 

For most of the J— and H— band observations, a 75 
line mm -1 grating in second order was used. This gave a re- 
solving power of 383 and 0.2 /im spectral coverage. For one 
case in the J— band, the 150 line mm -1 grating was used 
in third order, giving a resolving power of 1000 and a wave- 
length coverage 0.1 [im. For observations in the K— band the 
75 mm -1 grating was used in first order, giving a resolving 
power of 337 and 0.4 /jm spectral coverage. 

Standard nodding procedures were followed to subtract 
the sky background, with the slit oriented east-west on the 
sky. While observing, we nodded on to the detector so that 



we could integrate on the object while obtaining frames for 
sky subtraction. Bias frames were subtracted and the ob- 
servation images were then flat-fielded. An argon lamp and 
vacuum wavelengths are used here for the wavelength cali- 
bration. Correction for atmospheric extinction was applied 
by observing a spectral standard star at a similar airmass 
to the object. Flux calibration was achieved by extrapolat- 
ing J, H and K magnitudes from observed V magnitudes of 
flux standards according to known spectral type using the 
model of Koorneef (1983). 

We note that the measured counts in the "sky" frames 
were not always consistent with those in the "object" frames. 
From examination of the photometric consistency of the 
data we are led to conclude that any flux measurements 
have an uncertainty at the level of 10% (except for 3C 380 in 
which the uncertainty is ~20%). Data reduction was carried 
out using the CGS4DR package at the QMW and Edinburgh 
Starlink nodes. 

Figure 1 shows the observed spectra. The 14 SL sources 
are presented in RA order. The spectra of 3C 286 and the 
three serendipitous sources are given at the end. Note that 
there are only three SL sources where Ha was not detected. 
These are 3C 216, 4C 29.45, and 4C 56.27. In a further 
two sources, Ha is detected but is rather weak (3C 279 and 
B2 1308+32). In the remaining 9 a strong line was detected. 
Together with the 8 additional published sources that we 
use, this produces three sub-samples that we refer to : SL 
sources ( 22 objects), a sample of 19 SL sources with detected 
Ha emission, and a sample of 16 strong-EW sources. 



4.2 Emission line parameters 

Emission line measurements were performed using the 
DIPSO package on the ROE Starlink node. In Table 4 we 
present the observed Ha flux, rest frame EW and FWHM of 
the flux calibrated spectra. From these data both line and 
continuum luminosities (Lhc and L ) can be derived. The 
largest source of error comes from uncertainties in setting 
the local continuum (fitted by a first order polynomial). To 
minimize these errors each line was measured on at least 
three different occasions. In most cases the uncertainties are 
less than 10%. Experimentation with deblending the nar- 
row lines (Ha , [Nil], [SII] ) suggests that in most objects 
its presence does not introduce more than a 3% error into 
the flux of Ha . This maybe due to our inability to resolve 
and exclude the narrow lines from either sides of Ha . When 
the narrow Ha emission line is visible however, its contri- 
bution can change substantially FWHM. In these cases the 
measurements are corrected from narrow line contribution. 

In addition to our new data, Table 5 gives the 
Ha emission line parameters of 8 SL sources with pub- 
lished Ha spectra. These data are less homogeneous than 
the UKIRT spectra. We have searched for total line fluxes 
and EW (including also NLR) , but corrected FWHM (from 
NLR) values (references are given in Table 5). 
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5 COMMENTS ON INDIVIDUAL OBJECTS 

DA 193. The measured line flux and rest EW of the H/3 
emission line in this source are 1.7 1CP 14 erg cm -2 s _1 and 40 
A respectively. The total [OIII] flux is estimated as 9 1CP 15 
erg cm - s~ . The two spectra obtained of this source show 
an increase of the continuum flux at shorter wavelengths. 

4C 39.25. In the field of 4C 39.25 we have detected what 
appears to be an extended object (it occupies two detector 
rows) about 12 arcsec west of our source. We can identify 
this in the Schmidt plate as a point source which is brighter 
in the blue. From this and the spectrum (shown in Fig. 1) 
we conclude that it is a high redshift galaxy. 

3C 279. The broad feature at 1.1 /i of the spectrum is 
most likely due to the non-fully corrected absorption feature 
of the spectral standard. 

3C 334. The high resolution spectrum obtained in this 
source (see Fig 1) shows a structured Ha emission line pro- 
file but it has also a low signal-to-noise ratio. It is also 
likely not resolved in the continuum, and is arbritarily mul- 
tiplied in flux for comparison with the lower resolution 
Ha spectrum. 

3C 395. Two other objects, most likely stars, were de- 
tected in the slit at 3" east and at 15" west of the source. 
Their spectra are shown in Fig. 1. 



6 STATISTICAL RESULTS 

We have applied the Spearman rank correlation test (Press 
et al. 1988) to a set of nine parameters - to three emission 
line parameters (FWHM, EW, and Ln a ), to the optical 
luminosity at the Ha wavelength (L ), to three orientation 
indicators (R, /3 a and S ), to 9, and finally, to the Lorentz 
factor, 7. The resulting matrix of the correlation coefficient 
r s , and the two sided significance level of its deviation from 
zero, p, is shown in Table 6. The test was applied both to 
the sub-sample with detected Ha lines (N=19) and to the 
subsample with strong lines, divided at EW> 158 A (N=16). 

The strongest trend is seen between the Ha line and 
continuum luminosity (Lu a and L respectively), especially 
when only the strong-EW sources are considered. This is 
however a selection effect, the result of plotting two quanti- 
ties both multiplied by the same function of redshift. This 
trend is shown in Fig 2a where the three classes of objects 
CDP, CD, and LD are shown with separate symbols. (The 
three objects with Ha upper limits are also shown here). 
The strong- and weak-EW objects are very clearly distinct 
(see Fig 2b). The weak-EW objects have an EW at least 
a factor five less than the weakest strong-EW object. Note 
that the range of luminosities Lna is not very different in 
the strong- and weak-EW objects. As we shall see below, 
the FWHM does not show any such clear discontinuity, so 
that the two distinct trends in Fig 2 is almost certainly due 
to a difference in the continuum, not in the line, ft is also 



striking that all the weak-EW objects are polarised, ft seems 
then extremely likely that the weak-EW objects have a con- 
tinuum enhanced by relativistic beaming. Very likely there 
are two continuum components, the first of which is weak in 
the strong-EW objects, but is boosted to dominance in the 
weak-EW objects. The object 3C 345 with the largest EW 
amongst the CDP objects, looks like a good intermediate 
case, with significant but not dominant relativistic beam- 
ing. Broad band polarization and flux density spectra for 
this object do show the presence of two emission components 
(Wills 1991): one polarized, likely synchrotron jet radiation, 
and a second unpolarized possibly from a disc. Optical and 
Near-IR polarimetric and photometric monitoring also sup- 
port the two continuum emission components in this source 
(Smith et al. 1986). 

In all our subsequent plots, we reclassify objects as sim- 
ply strong-EW or weak-EW. The weak-EW class is identi- 
cal with the CDP class except for 3C 345, which has been 
classed as strong-EW. Within the strong-EW objects, we 
maintain a distinction between LD and CD objects. 

Equivalent width EW is correlated with the various ori- 
entation indicators, as indicated in Table 6 and Fig. 3, which 
collects these correlations graphically. (Note the inverted ab- 
scissa for the graphs versus 5 , R, and /3 a , chosen so that 
these can be easily compared with the correlation with 9.) 
Amongst the strong-EW objects, the strongest correlations 
are with S and with 9, in that order. The excellent corre- 
lation with S is as expected, as this indicates the amount 
of Doppler beaming. The fact that it is significantly bet- 
ter, with much less scatter, than the correlation with the 
indirect orientation indicator R , indicates that despite the 
uncertainties, stressed in Section 2, the derived values of S 
must be quite close to the correct values. The correlation 
with R is striking in a special respect as R is a kind of 
direct analogue of EW. While EW changes by only a factor 
of a few, R changes by a factor of several hundred. There 
is an anisotropy connected with the radio beaming, but it 
cannot be the same jet beaming effect (at least for the LD 
objects), but rather a second "weak beaming" effect. There 
is a very strong positive correlation between the line EW 
and the outflow angle of the radio jet in the whole sam- 
ple, but this is dominated by the cliff at small angles where 
the jet beaming sets in. The trend slightly weakens in the 
strong-EW sources and indicates grossly an increase of the 
EW by a factor of ~3 when the viewing angle turns from 5 
to 40 degrees. 

Next we look at the correlations of FWHM, which are 
collected in Fig 4. FWHM is correlated with all the orien- 
tation indicators, indicating a decrease of the emission line 
velocity field with increasing beaming. The direct correla- 
tion with 9 is suggestive of axisymmetric motions for the 
BLR. The most significant correlation seen is with viewing 
angle 9, as shown in Table 6. This is an indication that in 
deriving 9 we have done something meaningful. As further 
support for the reality of the 9 correlation, note that Table 
6 shows no correlation between FWHM and the jet Lorentz 
factor 7. The weak-EW objects fit well into the same trend 
as the strong-EW objects, being simply at one end of the 
trend, namely at small viewing angle. The trends weaken 
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somewhat when considering the strong-EW sample, which 
is however in a reduced range of 6. 



7 DISCUSSION 

Our main aim has been to compare Ha emission properties 
with observer viewing angle. The most exciting result is that 
both velocity and equivalent width do indeed correlate with 
orientation. The actual 6 derived from the SL motion plus 
Inverse Compton limit seems more reliable than the various 
indirect orientation indicators, as it gives stronger trends, 
and is physical. 

The FWHM vs 9 correlation is suggestive of a flattened 
structure for the line emitting material. Fig 5a tests this 
quantitatively, plotting V = Vo sin 9 with normalised line 
velocity of V = 13, 000 km s _1 at 0=90° . This is the disc fit 
of Wills & Browne 1986 to the H/3 in a larger AGN sample, 
extended to larger jet viewing angles. The figure also shows 
the effect of an isotropic component with V = 2, 000 km s _1 
(see for instance Brotherton et al. 1994, Corbin 1997) added 
in quadrature to the axisymmetric component. The agree- 
ment is good. The scatter is considerable, but we should cer- 
tainly of course expect some dispersion in Vo- How sensitive 
is this agreement with the predictions of an axisymmetric 
structure to the possible errors in 97 As described in section 
2, the value of 9 is strictly an upper limit, as the derived 
value of S is a lower limit. However, as argued there, S is 
very insensitive to the true Inverse Compton flux, and the 
correlations discussed in the previous section confirm that 8 
is unlikely to be in error by more than a factor of a few. We 
illustrate this sensitivity in Fig 6a by deriving the 9 values 
for each source with S increased by a factor three. It can be 
seen that the effect is still clear but the quantitative fit is 
not as good. It is also possible that the bulk speed is smaller 
than the pattern speed responsible for observed SL velocity. 
We do not illustrate the effect of this explicitly, but note 
that the net effect is similar to the S effect shown here. 

An important implication of a flattened structure for 
the broad line emitting region (BLR) concerns the mass M , 
of the black hole when this is determined from the line veloc- 
ity width assuming gravitational dynamics and isotropically- 
oriented motions (e.g. Wandel, Peterson & Malkan 1999, 
Kaspi et al. 2000). In a thin Keplerian BLR of inclina- 
tion angle equal to 9 the mass fit of a given line width is 
oc (sin 2 dVcosd)^ 1 which is ~ 1/ sin 2 9 at low 9 (Rokaki & 
Boisson, 1999). The inclination therefore may cause a sys- 
tematic underestimate of M , theoretically up to infinity. In 
more realistic orbital shape of the BLR, involving an aspect 
ratio of its flattening, the inclination may cause a system- 
atic underestimate of mass by a factor up to ~ 100 (Krolik 
2001). It is statistically implausible that the mass of the 
quasars is underestimated by such a large factor (this would 
happen if all sources are observed with 6—0) but it could be 
important for sources found accreting at super-Eddington 
rates (see Collin et al. 2002). The FWHM vs 6 correlation 
therefore, supports estimates of the black hole mass in which 



an orientation dependence of the line width is applied (e.g. 
Lacy et al. 2001, McLury & Dunlop 2002). 

The EW correlation shows the presence of two sepa- 
rate beaming effects - a strong effect presumably connected 
with relativistic beaming in the jet, and a second weaker 
effect, producing a factor ~ 3 change for 6 range from 5 to 
40 degrees. It is tempting to identify this second effect with 
the cos6 surface brightness effect expected from a flat ac- 
cretion disc. Fig 5b examines this possibility quantitatively. 
The dot-dashed curve shows the expected effect of Doppler 
enhancement for a jet with 7 = 30 (imposed mainly from the 
objects with large /? a ). This can explain the difference be- 
tween CDP objects and the rest, but not the general trend. 
The dotted line combines an isotropic continuum with a sec- 
ond separate jet-beamed continuum. This also does not fit 
the overall pattern. The dashed line shows the predicted ef- 
fect for a flat accretion disc, with standard limb darkening, 
where we should find EW = EW [1/3 cos 0(1 + 2cos0)] _1 . 
Fig 5b plots this with EW = 300 A. Note that apart from 
the normalisation there is no further degree of freedom. The 
solid line combines this disc-beamed continuum with a sec- 
ond separate jet-beamed continuum. The overall fit is good. 

In Fig 6b, we indicate the effect on 6 of changing 5 
by a factor of three, as explained in the FWHM discussion. 
The qualitative trend is still clear, but the quantitative fit 
is not as good. Given this uncertainty, we cannot hope, for 
example, to distinguish classes of accretion disc models such 
as flat discs, thick discs and flared discs. It is possible to 
restrict the range of 6 in the sample ( e.g. by decreasing H a ), 
more than in Fig 6b. In this case a jet model with a broad 
distribution of 7 could possibly account for the observed 
correlation, but requiring small 7 (<5) for lobe dominated 
objects (not seen in Table 2). 

Another worry is that we have implicitly assumed that 
the continuum light has a surface brightness effect with 6, 
but that the line emission does not. If the BLR is how- 
ever axisymmetric, as the FWHM effect seems to show, then 
Ha emission may come from the surface of the disc, and suf- 
fer a very similar effect with 6, removing the EW vs 6 effect. 
The effect would be the same as the dotted line in Fig 5b, 
showing isotropic continuum, which clearly does not fit the 
data. However, a likely scenario is that the line emitting 
material is in an axisymmetric structure but is not from the 
surface of an optically thick disc, being rather in a flattened 
system of clouds. Even if the lines are emitted from the sur- 
face of a disc, they will not necessarily suffer the same limb 
darkening law. 

If the observed optical continuum of the strong line class 
is emitted from a disc, it is not likely that it is markedly 
influenced by relativistic effects such as the gravitational 
redshift and lensing. These effects are stronger when the disc 
is viewed in the equatorial plane (and in the Kerr metric, 
Cunningham 1975), which is not the case for the SL objects. 
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8 CONCLUSIONS 

In this paper we have investigated the relation of the 
Ha emission line with radio beaming and outflow angle in a 
sample of SL sources. Most of the Ha spectra were observed 
by us at the UKIRT and are reported here. 

The sample is divided into two classes: 6 weak-EW and 
16 strong-EW quasars. The weak-EW class must have a rel- 
ativistically beamed continuum, since it has also high radio 
core dominance and high optical polarization. The strong- 
EW class shows a clear beaming effect, but one that cannot 
be relativistic beaming alone, since the line EW changes 
only by a factor of a few. The correlation with the actual 
viewing angle derived from the superluminal motions is very 
clear and nicely consistent with this weak beaming being due 
to the surface brightness effect expected from tilting a flat 
accretion disc. We also find that the FWHM of Ha is cor- 
related with the viewing angle in a manner quantitatively 
consistent with the idea that the line emitting material is an 
axisymmetric rotating structure. However, the line emission 
cannot come from the surface of an optically thick disc, oth- 
erwise the line beaming would cancel the continuum beam- 
ing. Overall it seems most likely that the optical continuum 
comes from a flat accretion disc whereas the line emission is 
from an axisymmetric system of clouds. 
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9 FIGURE CAPTIONS 

Figure 1 The observed flux calibrated spectra obtained 
with CGS4 at UKIRT . The 14 SL sources are presented 
in RA order. The spectra of 3C 286 and the three serendip- 
itous sources are given at the end. Note that there are only 
three SL sources where Ha was not detected (the expected 
wavelength is indicated in 3C 216, 4C 29.45, and 4C 56.27). 

Figure 2 (a): The Ha line versus continuum luminos- 
ity of the SL sample, (b): The EW versus the continuum 
luminosity. CDP, CD and LD objects are shown with cross, 
X and star symbols respectively. Small filled circle symbols 
show the UKIRT data. Upper limits are noted. 

Figure 3 The Ha EW versus R , /3 a , 5 , and 9, (in 
a,b,c and d respectively) of the 22 SL objects. In these 
and all subsequent figures, circles designate the strong-EW 
objects (open and filled symbols are core and lobe domi- 
nated sources respectively). Note the inverted abscissa for 
the graphs versus 8 , R, and /3 a , chosen so that these can be 
easily compared with the correlation with 9. 

Figure 4 The Ha FWHM versus the orientation in- 
dicators as in Fig 3, for the 19 objects with detected line 
emission. 

Figure 5 (a): FWHM versus 9 compared with models. 
The dashed curve shows an axisymmetric component (V = 
13, 000 sin 9 km s _1 ) and the solid curve gives its quadrature 
sum with an isotropic one (V=2,000 km s -1 , the dotted 
curve), (b): EW versus 9 compared with models. Note the 
inverted EW axis indicating the brightness of the underlying 
continuum. The dot-dashed curve shows the expected effect 
of Doppler enhancement for a jet with 7=30. The dashed 
curve shows the predicted effect for a flat accretion disc with 
a standard limb darkening, and the solid curve combines 
disc and beamed component. The dotted curve combines an 
isotropic and a beamed component. 

Figure 6 Same as Fig 5, but with 9 recalculated using 
8 bigger by factor three to show sensitivity to errors. 
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Figure 1. 
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Figure 1. - continued 
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Table 1. Names, redshift, SL speed and radio core dominance of the SL sample taken from VC94. The observed radio core flux and 
angular size at a given wavelength, and the X-ray flux at 1 keV are given together with the references. 
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Class 


z 
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S„ 
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o d 
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(inn7 i 1 na 
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TTT 
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1.25 
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A A 
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A Q 

43 
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T3 
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0133+20/ 




A T 

4/ 


T 





425 


a n 
4.9 


— 1 


30 


0. 


084 
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0.15 


V93 
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A A 

44 


rJ 


U4oU+U0z 


on 
QKj 


1 on 




U. 
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K A 
0.4 
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.Oil 


o 
O. 


linn 


o.u 


n Af\ 
U.4U 


T5Q 1 

rol 


1U. 


UU 


rSlvl 


0552+398 




193 


CD 


2 


365 


2.3 





an 


z. 


620 


a a 
8.4 


n 7q 
0.73 




n 
U. 


An 
.49 


T > 


0850+581 


4C 


58.17 


CD 


1 


.322 


5.2 





.27 


0. 


940 


5.0 


0.48 


PR 





.06 


B 


0906+430 


3C 


216 


CDP 





669 


5.1 


-0 


.01 





880 


5.0 


0.10 


PR 





11 


B 


0923+392 


4C 


39.25 


CD 





.699 


5.3 


1. 


.30 


0. 


200 


22.2 


0.05 


A97 





.58 


B 


1040+123 


3C 


245 


CD 


1 


.029 


4.1 





.00 





590 


10.7 


0.33 


HR 





.13 
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1156+295 


4C 


29.45 


CDP 





729 


34.8 





.83 


1. 


400 
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0.12 


H90 
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1226+023 


3C 


273 


CD 
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10.7 
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1. 


820 
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0.21 


MOO 


10 


.89 
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1253-055 


3C 


279 


CDP 





538 


12.2 


1 


.10 


12. 


310 


22.0 


0.08 


W01 





.96 


B 


1308+326 


B2 




CDP 





.996 


27.7 


1 


.60 


0. 


540 


5.0 


0.18 


G93 





.13 


B 


1618+177 


3C 


334 


LD 





555 


2.5 


-0 


.64 





.086 


10.7 


0.20 


H92 





.20 
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1641+399 


3C 


345 


CDP 





595 


12.4 
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22.0 


0.13 


RO 





.57 
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4C 


34.47 


LD 
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-0 
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0.24 


HI 
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.37 


B 


1823+568 


4C 


56.27 


CDP 
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3.4 
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.30 


0. 
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0.31 
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B 


1830+285 
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H2 





.31 
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1. 


600 


5.0 


0.35 


M87 
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1. Unless otherwise noted /3 a is from VC94 

2. /3 a and R from BOO and WB 

3. /3 a and R from A96 and WB 



REFERENCES: 
A96: Alef et al. 1996 
A97: Alberdi et al. 1997 
B: Brinkmann et al. 1997 
BOO: Brunthaler et al. 2000 
B95: Brinkmann et al. 1995 
BM: Bloom & Marscher 1991 
C90: Chariot 1990 
G93: Gabuzda et al. 1993 



G94: Gabuzda et al. 1994 
HI: Hooimcycr et al. 1992 
H2: Hooimeyer et al. 1992 
H90: McHardy et al. 1990 
H92: Hough et al. 1992 
HR: Hough & Rcadhcad 1987 
L99: Lara et al. 1999 
LO: Ledden & Odell 1985 



MOO: Mantovani et al. 2000 
M87: Madau et al. 1987 
P81: Pauliny-Toth et al. 1981 
PR: Pearson & Readhead 1988 
PW: Polatidis & Wilkinson 1998 
RO: Ros, Zcnsus & Lobanov 2000 
V93: Vcrmculen et al. 1993 
W01: Wehrle et al. 2001 
WB: Wills & Browne 1986 
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Table 2. The derived <5 , 7 and 6 (in degrees) of the SL sources listed in increasing 9. 



Name 




s 




7 




(9° 


3C 279 


73. 
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75 





.25 


3C 216 


42 


.58 
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.61 





.32 


4C 29.45 
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.22 
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.55 
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1308+326 


11 


59 
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.52 


3C 395 
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.26 
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.00 
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.45 
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.33 
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3C 380 
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.22 


13 
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,72 
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.99 
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.31 
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11 


3C 47 
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25. 


37 


22 


.75 


3C 245 


1 


.73 


6 


.12 


23. 


.41 


4C 28.45 





.83 


8. 


.02 


31. 


.01 


4C 34.47 





.19 


29 


24 


35 


05 


3C 334 





.40 


9 


.19 


42 


68 


3C 390.3 





.45 


7 


85 


43 


.76 



Table 3. The UKIRT observing log. Unless otherwise noted, the observations were made in J band. 



Name 


Obs. date 


Seeing 
arcsec 


Exp. time 
min 


Notes 


DA 193 


29/3/94 


1.16 


46 


K, grating in first order 




29/3/94 


1.35 


72 


H, 600 km s" 1 


4C 58.17 


29/3/94 


1.32 


55 


H, 600 km s- 1 


3C 216 


29/3/94 


1.31 


27 




4C 39.25 


27/3/94 


1.07 


60 


2 objects in the slit 


3C 245 


28/3/94 


1.10 


80 




4C 29.45 


28/3/94 


1.04 


35 




3C 279 


27/3/94 


1.23 


88 




1308+326 


27/3/94 


1.02 


86 




3C 286 


29/3/94 


1.06 


75 


z=0.849 


3C 334 


27/3/94 


1.04 


60 






28/3/94 


1.30 


85 


150 1/mm, R=1000 


3C 345 


27/3/94 


1.08 


60 




4C 56.27 


29/3/94 


1.38 


30 




3C 380 


28/3/94 


1.27 


40 




4C 28.45 


28/3/94 


1.06 


15 




3C 395 


29/3/94 


1.15 


30 


3 objects in the slit 
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Table 4. The observed Ha flux (in lCT 14 erg cm- 2 s" 1 ), rest frame EW and FWHM of the UKIRT Ha spectra. 



Name 


Flux 


EW 


t WHM 






A 


Km s 


DA 193 


6.00 


240 


2450 


4ty oo.l / 


q An 


Z4U 


a 1 nn 
41UU 


3C 216 


<0.05 


<10 




oy.zo 


q fin 


one 


4UUU 


3C 245 


0.84 


316 


3250 


4C 29.45 


<0.06 


<8 




3C 279 


1.40 


8 


1400 


1308+326 


0.90 


30 


2350 


3C 286 


4.6 


254 


2600 


3C 334 


13.60 


446 


6150 


3C 345 


3.25 


186 


3600 


4C 56.27 


<1.33 


<2 




3C 380 


3.37 


380 


3500 


4C 28.45 


6.67 


354 


4050 


3C 395 


28.60 


158 


2050 



Table 5. The observed Ho flux (in 10 14 erg cm 2 s 1 ), rest frame EW and FWHM and the references to the published Ha spectra. 



Name Flux EW FWHM Ref. 

A Kins" 1 



III Zw 2 


136 


445 


3600 


N79.R85 


3C 47 


10.5 


392 


8400 


JB 


3C 120 


77.0 


299 


1846 


RB 


3C 273 


650.0 


369 


3260 


JB 


4C 34.47 


65.0 


589 


2300 


S81, EH 


3C 390.3 


80.6 


519 


11900 


N82,L96,EH 


4C 73.18 


38.0 


365 


3100 


JB,EH 


BL Lac 


28.7 


7.3 


4000 


C96 



REFERENCES: 

C96: Corbett et al. 1996, EH: Eraclcous & Halpcrn 1994 
JB: Jackson & Browne 1991a, measured in scanned spectra 
L96: Lawrence et al. 1996, N82: Netzcr 1982 
N79: Ncugebauer et al. 1979, R85: Rafanelli 1985 
RB: Rokaki & Boisson 1999, S81: Soifer et al. 1981 
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Table 6. The matrix of the correlation coefficient r s , and the two sided significance level of its deviation from zero, p, for the sub-sample 
with detected Ho lines (N=19) and for the subsample with strong lines (N=16). Correlations with p < 0.02 are shown in bold type. 



FWHM EW L Ln. 





r S 


P 


rs 


V 


rs 


P 


r S 


P 






Strong- and weak-EW objects, N= 


— 1 o 






EW 


0.37 


0.115 














L 


-0.28 


0.244 


-0.42 


0.074 












0.02 


0.949 


0.10 


0.687 


0.71 


0.001 






R 


-0.38 


0.107 


-0.62 


0.005 


0.08 


0.744 


-0.29 


0.228 


/3a 


-0.46 


0.047 


-0.59 


0.008 


0.37 


0.115 


-0.01 


0.983 


5 


-0.47 


0.043 


-0.78 


7.5 10" 5 


0.28 


0.241 


-0.22 


0.363 


1 


-0.32 


0.178 


-0.02 


0.930 


0.42 


0.075 


0.17 


0.477 


Q 


0.57 


0.011 


0.70 


0.001 


-0.49 


0.033 


-0.01 


0.955 








Strong-EW objects, 


N=16 








EW 


0.33 


0.214 














L 


-0.07 


0.805 


-0.52 


0.040 












-0.14 


0.615 


-0.34 


0.192 


0.96 


5 10- 9 






Ft 


-0.41 


0.113 


-0.44 


0.092 


0.09 


0.747 


0.01 


0.976 


/3a 


-0.33 


0.215 


-0.52 


0.037 


0.25 


0.345 


0.15 


0.572 


8 


-0.33 


0.207 


-0.74 


0.001 


0.17 


0.535 


0.01 


0.983 


7 


-0.11 


0.674 


0.13 


0.623 


0.22 


0.412 


0.26 


0.333 


e 


0.46 


0.076 


0.67 


0.005 


-0.40 


0.122 


-0.28 


0.289 
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Figure 2. 
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Figure 3. 
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Figure 4. 
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A Disc in the SL Quasars? 
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Figure 6. 
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